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ABSTRACT 


This  report  describes  the  research  accomplished  in  the  second 
quarter  of  the  program  entitled  "CO  Chemical  Laser  Research". 

Briefly,  the  objectives  of  this  research  are  to  experimentally  study 
single  line  operation  of  a CO  chemical  laser,  and  to  develop  a 
theory  of  cw  CO  chemical  lasers. 

The  first  task  of  the  program,  a demonstration  of  single  line 
operation  of  a cw  CO  chemical  laser,  was  completed  in  the  first 
quarter  and  reported  in  the  First  Quarterly  Progress  Report 
(23  February  1972  - 1 June  1972). 

The  theoretical  treatment  of  the  cw  CO  chemical  laser  (Task  3) 
is  the  main  subject  of  this  report.  In  order  to  investigate  the 
collisional  relaxation  effects  occurring  in  CO  chemical  laser  media, 
the  master  relaxation  equations  for  CQ-M  systems  were  solved  by 
computer  calculations.  The  relaxation  of  CO*  from  the  0 + CS  chemical 
reaction  distribution,  was  studied  for  the  cases  of  CO*,  CO*-He, 
C0*-02,  CO*-CO ( 300  K),  and  C0*-02  + CO(300  K) . Unless  a large 
excess  of  some  collision  partner  is  present  (CO*:M  of  1:100  or 
more),  the  dominant  collisional  relaxation  phenomena  will  be  CO*-CO* 
or  C0*-C0(300  K)  processes. 

The  complete  theory,  including  fluid  flow,  chemical  kinetics, 
relaxation,  and  optical  properties,  is  described  in  detail.  The 
zero  optical  intensity  theory  is  used  to  model  the  MDRL  CO  one ui cal 
laser  for  two  cases  of  interest:  with  and  without  CO (300  K/  injec- 
tion into  the  reacting  stream.  The  results  agree  in  general  with 
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experiment,  but  predictions  of  R-branch  gains  from  the  theory  suggest 
that  the  published  CO*  distribution  formed  by  the  0 + CS  pumping 
reaction  has  a limited  range  of  validity. 
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I . INTRODUCTION 


This  report  describes  progress  made  in  the  second  quarter 
(1  June  1972  - 1 September  1972)  on  the  program  "CO  Chemical  Laser 
Research"  at  the  McDonnell  Douglas  Research  Laboratories  (MDRL) . 

The  technical  objectives  of  this  contract  are  as  follows: 

Task  1.  Experimental  demonstration  of  single  line 

operation  with  a cw  CO  chemical  laser  based  on 
the  0 + CS2  reaction. 

This  task  was  completed  in  the  first  quarter  and  the  results  reported 
in  the  First  Quarterly  Contract  Report  (23  February  1972  - 
1 June  1972) .1 

Task  2.  Experimental  study  of  the  effects  of  various 
parameters  on  single  line  output  power  for 
the  chemical  laser  of  Task  1. 

This  task  is  in  progress  and  results  will  be  reviewed  in  subsequent 
reports . 

Task  3.  Formulation  of  a theory  of  the  cw  CO  chemical 

laser  to  calculate  both  zero  intensity  parameters 
and  the  interaction  of  the  stimulated  optical 
fields  with  the  active  medium. 

The  theory  of  cw  CO  chemical  lasers  has  been  completed  for  the  zero 
intensity  case.  This  theory  and  results  obtained  from  computer 
calculations  with  it  will  be  the  main  subject  of  the  present  report. 

The  principal  investigator  is  W.  Q.  Jeffers.  R.  F.  Webbink 
and  F.  S.  Skinner  also  participated  in  this  effort  in  the  second 
quarter. 
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II.  PROGRESS  AND  RESULTS 

Task  3 

2.1  Collisional  Relaxation  of  CO* 

Our  effort  on  the  theoretical  description  of  cw  CO  chemical 
lasers  was  directed  to  two  major  problems  in  the  second  quarter.  In 
the  first  problem,  we  have  programmed  the  master  relaxation  equations 
for  CO-M  mixtures.  The  program,  called  RELAX  1,  takes  an  arbitrary 
initial  distribution  of  CO  in  excited  states  up  to  v = 50  and  cal- 
culates the  subsequent  distributions  which  evolve  by  single  quantum 
V-V,  V-T,  and  spontaneous  emission  processes.  The  V-V  and  V T 
processes  are  included  for  CO-CO  and  CO-M  collisions  where  M is 
either  a monatomic  or  diatomic  species.  In  the  case  where  M is 
diatomic,  the  M distribution  is  assumed  to  be  in  thermal  equilibrium 
at  the  gas  kinetic  temperature  which  is  constant  at  300  K. 

Transition  probabilities  for  CO-M  collisions  (M  = CO,  He,  02) 
for  both  V-V  and  V-T  energy  transfers  have  been  calculated  by  our 
previously-developed  techniques.2'3  The  V-V  calculations  for  CO-CO 
exchanges , 

C0(v)  + CO(v')  ■*  CO(v  ± 1)  + CO  ( v ' + 1)  (-1) 

have  been  done  for  both  loig-range  and  short-range  interactions, 
and  for  (v,v')  < 50.  The  transition  probabilities  from  the  short- 
range  interactions  are  limited  to  values  of  0.125  by  the  theory.3 
The  long-rang  interactions  are  treated  by  a perturbation  theory 
result,  and  we  have  set  an  upper  limit  on  the  calculated  probabil- 
ities, P,  by  using  a pseudo-probability  P'  defined  by 
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p.  = I (2) 

1 f eP  (e  = 2.71828  ...) 

The  V-T  exchange  with  He  was  treated  using  the  standard  SSH 
theory4  of  V-T  transfer.  The  calculated  results  for 

CO{v)  + He  -*•  CO(v  - 1)  + He  (3) 

were  adjusted  to  experimental  data5  by  a scale  factor.  The  exchange 
with  02  was  also  treated  in  the  relaxation  equations  as  a V-T  col- 
lision, but  the  probability  for  this  event, 

CO  (v)  + 02(0)  CO  ( v - 1)  + 02(1) 

3 

was  calculated  from  the  V-V  theory  for  short-range  interactions. 

5 

Again,  the  computed  results  were  adjusted  to  experimental  data  by 
a scale  factor. 

The  first  calculation  done  was  the  case  for  chemically-formed 

6 

CO  relaxing  from  the  initial  distribution  of  the  0 + CS  reaction. 
This  reaction  is  well  established  as  the  pumping  step  in  CO  chemical 
lasers, ^ and  its  distribution  has  been  deduced  from  at  least  two 
different  experiments.5'6  The  relaxation  calculation  includes  CO-CO 

• • 8 —i  « i 

V-V  and  V-T  processes  as  well  as  spontaneous  emission.  inis  cal- 
culation assumes  that  n,  the  total  density  of  CO,  is  3.215  x 1015  cm 
and  is  distributed  in  the  excited  states  as  shown  in  the  t = 0 curve 
of  Fig.  1.  This  distribution  represents  a slight  smoothing  of  the 
distribution  obtained  experimentally  by  Hancock  and  Smith.  Figures 
1-5  display  the  results  as  functions  of  both  vibrational  quantum 
number  v and  time  t. 
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Fig.  1 Relaxation  of  chemically- formed  CO*  at  P = 0.1  Torr,  T - 300 
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In  general,  relaxation  from  the  peaked  chemical  distribution 
is  rapid.  Population  spreads  to  both  higher  and  lower  levels  and 
tends  to  form  a flat  distribution  at  intermediate  times  (t  £ 1000  psec) . 
At  longer  times  the  distribution  could  be  described  by  a low  vibra- 
tional temperature  at  low  v (v  ~ 6)  and  a high  vibrational  temperature 
at  high  v (v  > 7) . Figure  2 is  a plot  of  maximum  (with  respect  to 
J)  P-branch  gains  ap  which  exist  on  each  vibrational  transition. 

There  is  a rapid  decrease  in  the  gains  on  the  bands  directly  excited 
by  the  chemical  reaction.  Bands  at  low  and  high  vibrational  levels 
build  up  gain  as  the  population  distribution  spreads,  and  subsequently 
relax  back  toward  smaller  gains.  R-branch  gains  ctR  exist  only  on  the 
positive  slope  side  of  the  chemical  distribution  (Fig.  3) , spread 
toward  low  v,  and  decrease  rapidly.  Even  though  V-V  transfer  removes 
no  quanta  from  the  system,  the  redistribution  decreases  the  optical 
gain  coefficients.  All  of  these  effects  can  be  seen  again  in  Figs.  4 
and  5 in  the  time-dependence  of  the  maximum  P-  and  R-branch  gains.  It 
is  clear  that  CO-CO  V-V  transfer  is  a rapid  mode  of  redistribution  of 
population  from  the  chemically-formed  distribution. 

In  order  to  evaluate  the  effects  of  various  collision  partners 
which  would  be  present  in  a real  CO  chemical  laser,  calculations 
were  run  for  the  chemical  distribution  again  with  n£0|.a]_  = 3.215  x 
10 1 5 cm- 3 (corresponding  to  P - 0.1  Torr)  and  with  1.0  Torr  of  He,  O2, 
and  "cold"  CO.  The  calculations  were  run  independently  to  separate  the 
effects  of  the  various  collision  partners.  Figures  6 and  7 show  the 
results  of  the  CO*: He  (1:10)  calculation.  Both  the  distributions 
and  gain  coefficients  are  negligibly  different  from  the  case  with 
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Fig.  5 Time  dependence  of  R-branch  maximum  gain  transitions  for  v - 5,  10. 


Relaxation  of  CO*:IIe  (1:10)  at  PCQ*  = 0.1  Torr,  PRe  = 1.0  Torr 
T = 300  K. 


Fig.  7 Maximum  P-branch  gain  coefficients  ap  resulting  from  the  time- 
dependent  distributions  of  Fig.  6. 


CO*  (CO  initially  in  the  chemical  distribution)  alone.  Thus,  we 
conclude  that  the  only  effect  of  He  in  CO  chemical  lasers  is  to  act 
as  a temperature-moderating  diluent  (see  below). 

Figures  8-10  are  the  results  of  the  C0*:02  (1:10)  calculation. 

The  energy  separation  of  the  v = 1 to  v = 0 states  of  02  is  1556.2  cm  , 
which  will  be  nearly  resonant  with  Av  = 1 CO  transitions  at  v = 24. 

The  V-V  transfer  probability  for  the  collisions  of  Eq . (4)  thus  goes 

through  a maximum  with  respect  to  v at  v = 24.  Therefore,  a large 
excess  of  0?  with  CO*  will  act  as  a barrier  to  V-V  pumping  of  CO 
population  beyond  v = 24.  Figure  8 shows  this  effect  quite  clearly 
at  times  of  1000  ysec  or  longer,  where  with  02  the  populations  in 
levels  v > 21  are  progressively  lower  than  the  case  with  only  CO*. 

The  CO  P-branch  gain  coefficients  for  v £ 20  are  larger  by  typically 
35%  at  relatively  long  times  ('v-  3000  jisec  at  0.1  Torr  CO*,  1.0  Torr 
02)  because  of  suppression  of  upward  V-V  pumping  by  02 . Therefore, 
with  C0*:02  ratios  of  1:10  or  more,  the  collisional  effects  of  02  in 
the  CO  chemical  laser  prevent  upward  V-V  pumping  of  CO  beyond  levels 
of  about  v = 20. 

Figures  11-15  concern  the  case  of  "cold"  CO  and  CO*.  We  have 

observed  in  CO  chemical  lasers  that  the  addition  of  vibrationally 

9 

cold  CO  increases  the  output  power  typically  by  a factor  of  2.  In 

CO  pulsed  electric  discharge-  lasers,  the  differential  relaxation 

10 

effect  of  "cold"  CO  is  the  mechanism  which  produces  inversion. 

As  shown  in  Fig.  11,  this  differential  relaxation  of  the  CO*  dis- 
tribution tends  to  increase  the  inversion  ratios  for  v > 4.  At  the 
same  time  the  magnitudes  of  all  of  the  excited  populations  decrease 
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Fig.  15  Time  dependence  of  R-braneh  gain  coefficients  from  Fig.  11. 


more  rapidly  due  to  the  predominance  of  the  de-exciting  collisions, 


CO  (v)  + CO(O)  -*  CO  ( v - 1)  + CO(l).  (5) 

The  result  is  a dip  in  the  population  distribution  which  decreases 
in  magnitude  and  propagates  toward  high  v as  the  relaxation  pro- 
gresses. In  Fig.  11  the  minimum  of  this  dip  an  be  seen  at  v - 4 
at  100  Msec,  v = 5 at  300  usee,  v = 8 at  1000  Msec,  etc.  At  some 
given  time  in  the  evolution  of  the  CO*  distribution,  there  will  be 
a range  of  v for  which  the  gain  coefficients  are  larger  with  "cold" 

CO  than  without  it.  For  example  at  100  Msec,  the  ap  are  larger  at 
v = 7,  8,  9,  10  in  Fig.  12  than  in  Fig.  2.  For  v > 10  the  ap  are 
nearly  the  same  for  the  cases  CO*  and  CO*: CO  (1:10),  and  at  v < 7 
the  "cold"  CO  decreases  ap.  The  effects  on  the  R-oranch  gains  can 
be  seen  most  clearly  by  comparing  Figs.  5 and  15.  With  CO*  alone 
the  chemical  distribution  predicts  large  total  inversion  ratios  at 
t = 0 for  the  bands  from  12  -*•  11  downwards.  As  the  CO*  distribution 
spreads  to  lower  v by  relaxation,  the  high  bands  (e.g.,  10  -*  9)  show 
a monotonic  decrease  in  a . Lower  bands  (e.g.,  5 4 ) rise  to  a maxi- 

K 

mum  and  then  decrease,  since  the  magnitudes  of  the  initial  populations 
are  small.  With  the  addition  of  "cold"  CO,  the  differential  relax- 
ation tends  to  inhibit  relaxation  away  from  total  inversion.  This 
can  be  seen  in  Fig.  15  where  the  10  -*■  9 aR  comes  back  to  a second 
peak  at  500  Msec. 

p 

Figure  16  shows  plots  of  the  decay  of  a10  (the  maximum  P-branch 
gain  on  the  v = 10  -*■  v = 9 transition)  for  the  cases  mentioned.  The 
choice  of  coordinates  in  this  figure  is  governed  by  the  linear 
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dependence  of  gains  on  the  density  of  excited-state  CO  in  the  region 
where  pressure-broadening  is  negligible,  and  by  the  linear  dependence 
of  the  logarithmic  derivatives  of  the  level  populations  on  gas  density 
in  regions  where  collisional  processes  dominate  the  excited  state 
population  evolution.  The  time  evolution  of  a?0  may  then  be  deter- 
mined from  Fig.  16  over  a reasonably  wide  range  of  pressures: 

10-2  Torr  < PCQ*  < 1 Torr.  For  tP  < 0.5  ysec-Torr,  all  curves 

converge  on  otfo/P^Q*  = 3.7  cm  JTorr 

In  summary,  we  have  now  examined  the  collisional  relaxation 
effects  of  those  species  which  will  be  of  major  importance  in  0/CS2 
CO  chemical  lasers.  The  predominant  effects  will  certainly  be 
either  CO*-CO*  relaxation,  or  CO*-"cold"  CO  relaxation.  Such  col- 
lision partners  as  He  and  O2  will  have  only  slight  effects  on  the 
chemically-formed  CO. 
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2.2  Concurrent  Fluid  Flow,  Chemical  Kinetics,  and  Collisional 
Relaxation  of  CO* 

The  program  to  describe  the  combined  effects  of  flow,  chemistry, 
and  relaxation  has  been  set  up  to  model  the  flow  and  geometry  of  the 
MDRL  chemical  laser. The  chemistry  and  relaxation  have  been  treated 
in  a completely  general  way,  and  could  be  adapted  to  any  particular 
fluid  flow  situation. 

Figure  17  shows  a cross  section  of  the  reactant  flow  tubes  used 
in  the  experimental  laser  device.  Some  initial  set  of  species  [A] 
is  contained  in  the  main  (round)  injector  at  pressure  PA  and  temper- 
ature T . This  gas  is  pumped  through  a slot  of  width  and  length 
L into  the  cavity  region  which  is  assumed  to  be  at  pressure  Pcavity* 

No  pressure  gradients  are  allowed  within  the  cavity  region.  The 
initial  stream  is  assumed  to  be  non-reactive.  The  expansion  from  the 
injector  interior  to  the  cavity  region  is  assumed  to  be  ieen tropic , 
and  its  temperature,  thickness,  and  flow  velocity  are  computed  using 
this  assumption.  A specification  of  the  flow  rates  F^  of  all  of  the 
species  in  [A]  is  given  and  all  species  are  assumed  to  behave  as 
ideal  gases.  For  x < -xfi  the  initial  stream  quantities  are  given  by 


n 


v 


Pcavity 

kT 


cavity 

PA 


) 


1/2 


(7) 


(8) 
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2 — / 
nvL 


(9) 


where 

g^  -•  1/2  times  the  average  number  of  degrees  of  freedom  per 
particle  in  the  mixture  [A] , excluding  vibration.  This 
represents  the  heat  capacity  (excluding  vibration)  of 
an  ideal  gas  at  constant  volume, 
m^  = average  mass  per  narticle  in  [A] 
n = particle  number  density 

v = flow  velocity  of  the  initial  stream  in  the  x direction 
z = thickness  of  the  initial  stream  (of  length  L) 

F = total  fl:w  rate  of  the  initial  stream  (particles/sec). 

/a 

All  of  the  stream  quantities  v,  z,  and  T are  constant  until  the 
initial  stream  reaches  the  first  injector.  This  follows  from  the 
assumptions : 

1.  The  initial  stream  is  non-reactive,  hence  no  changes  in 
internal  energy  or  chemical  species. 

2.  The  gas  is  assumed  inviscid  and  does  not  interact  with 

the  surrounding  gas  other  than  to  adjust  its  internal 

pressure  to  p . . . 

c ^cavity 

3.  There  are  no  pressure  gradients  in  the  cavity,  hence 

32-0. 

dx 

At  x = -x  the  initial  stream  reaches  the  first  transverse 
injector  pair  B.  The  flow  of  reactants  [B]  from  the  injector  pair 
is  treated  in  the  same  manner  as  the  expansion  from  the  main  injector 
A,  with  the  use  of  the  parameters  [B] , p0,  Tg , wfi,  L,  an^  Fg.  The 
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flow  from  the  injector  pair  B is  assumed  to  be  in  the  z-direction, 
so  the  injectors  are  properly  "transverse"  injectors.  At  x = -x_, 
the  stream  quantities  are  recomputed  under  the  assumption  that  the 
combined  flows  [A]  and  [B]  mix  instantly  without  reaction  and  come 
to  an  internal  pressure  pc. 

In  order  to  compute  the  stream  quantities  at  x = -xA,  the  fol- 
lowing conservation  laws  are  used: 

Conservation  of  particles: 

F = Fa  + F (10) 

Conservation  of  mass: 

n;  = faSa  + fbSb  • (11) 

Conservation  of  number  of  degrees  of  freedom: 

Pg  = FAgA  + Vb  • (12) 


Conservation  of  x-momentum: 

Fmv|  = FAmAv  j , (13) 

X=-XB  X<-Xfi 

Conservation  of  energy  (at  constant  pressure) : 

(fa+  fb)  (ySv2  + gkT  + pi)  | 

x=-xB 

= FA(TV2  + ^AkT  + pn)  I (14) 

X<-xb 


x Fn(4Snv2n  + gRkTR  + 


B 2 B B 


B B 


The  quantities  Vg,  gB,  Tg,  and  ng  are  evaluated  for  the 
expanded  stream  from  injector  B be  fore  mixing  with  the  main  stream. 
vn  is  in  the  z-direction.  Thus,  at  x = -xR  the  gases  have  "mixed" 
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(for  fluid  flow  purposes)  and  the  new  chemical  species  set  [A  + B]  is 
known,  along  with  the  new  values  g,  m,  T,  n,  v,  and  z.  At  x - -x^ 
the  chemical  reactions  begin,  but  no  instantaneous  changes  can 
occur.  To  model  the  real  situation  of  finite  mixing  times,  the 
concentrations  of  one  reaction  set,  [B] , are  made  available  to  the 
chemical  reactions  with  [A]  as 


dYi 

dx 


n . -y . 
1 

o . 

^mix 


^i  dz 
z dx 


(15) 


i. 

where  ni  is  the  density  of  the  itn  species  in  [B] , and  yi  is  that 
portion  of  n.  which  has  been  made  available  to  reactant  set  [A] . 
The  second  term  arises  from  stream  expansion,  z being  the  stream 


thickness.  The  mixing  length,  £ ix,  must  be  furnished  as  an  input 
parameter.  The  quantity  £mix  might  be  calculated  on  the  basis  of 
some  model  of  the  flow  and  mixing  situation  (such  as  diffusive 


mixing) , or  it  might  be  chosen  as  an  experimental  estimate. 

The  chemical  kinetics  have  been  limited  to  three  reactions 

since  collisional  effects  on  CO  are  to  be  allowed  with  all  reactant 

and  product  species.  In  addition,  there  is  evidence  that  the 

reactions  below  are  a fairly  complete  description  of  O/O2/CS2 

7 

chemistry  at  temperatures  below  500  K: 

k 


• 1 6 


0 + CS2  -*•  C.S  + SO  - 30  kcal/mol 

k16  « (1.07  x ic-10,-664.67  hc/kT  12 


(16) 
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(17) 


and 


1 7 


0 + CS  CO*  + S - 75  kcal/mol 

k17  = (4.20  x 10-io}  “709.01  hc/kT , 


13,14 


k i 8 

S + O2  -*•  SO  + 0-  6 kcal/mol 

, /->  -0.0  hc/kT. 

k ! 8 = (2 . 00  x 10  ) e 


15 


(18) 


The  rates  given  above  are  for  unit  species  densities,  and  the  acti- 
vation energies  are  in  units  of  cm  * molecule  1 . References  refer 
to  the  sources  for  the  rate  constants;  for  the  0 + CS  reaction,  an 
activation  energy  was  deduced  from  a fit  to  rate  constants  at  1100  K 
and  300  K.13 

To  calculate  stream  properties  at  x > -x_  the  chemical  kinetics 
must  be  included  since  the  chemical  species  are  changing  and  the 
reactions  release  heat.  As  before,  the  stream  is  allowed  to  expand 
or  contract  so  that  its  internal  pressure  is  Pcavity‘  Let  the 
&th  reaction  be  described  by  a rate  k^  per  unit  mass,  energy 

release  e0,  and  r.  the  number  of  particles  of  species  i created 
or  destroyed  in  each  elementary  step  of  the  reaction.  Then  the 
differential  equation  for  the  temperature  of  the  stream  is 


where 


(19) 


^rad 


net  energy  loss  through  radiative  processes  per 
unit  mass 
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U = energy  stored  in  vibration  per  unit  mass  (only  CO) 
k = Boltzmann  constant. 

The  stream  thickness  z is  calculated  from  the  differential 
equation 

dz  _ — \y enk0dx  - vdQ  . - vdU 

a v(g+l)kT[j-^  A A rad  v 

(20) 

- - g>2>uk<;dx]  • 

i i J 

Between  the  first  transverse  injector  pair  B and  the  second  trans- 
verse injector  pair  C,  the  stream  quantities  T and  z are  calculated 
from  Eqs . (19)  and  (20),  the  flow  velocity  is  constant  after  the 
discontinuity  at  x = -x  , and  the  chemical  species  are  calculated 
from  the  rate  equations  described  by  the  reaction  steps  (16),  (17), 

and  (18) . 

After  the  first  injection  point  x = -xB,  reaction  17  begins  to 
form  CO*.  The  initial  populations  of  excited  states  are  usually 
taken  to  be  zero,  however,  the  program  also  allows  for  an  initial 
Boltzmann  distribution  at  temperature  T.  The  probability  of  for- 
mation by  reaction  17  into  each  excited  state  v is  the  normalized 
initial  distribution  of  Fig.  I.6  As  in  the  relaxation  calculations 
the  distribution  of  CO*  is  calculated  by  integrating  the  master 
equations.  Both  V-T  and  V-V  collisional  energy  exchange  processes 
are  allowed  for  each  chemical  species  existing  in  the  set  of 
reactions  16,  17,  and  18. 
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For  the  V-V  processes  the  following  collisions  are  considered: 


CO  (v)  + CO(v')  - CO(v  ± 1)  + CO  ( v ' + 1)  (21) 

CO  (v)  + CS2(00°0)  ->  CO(v  - 1)  + CS2(00°1)  (22) 

CO  (v)  + CS  (v  = 0)  CO  (v  - 1)  + CS  (v  = 1)  (23) 

CO(v)  + SO  (v  = 0)  -►  CO  (v  - 1)  + SO  ( v = 1)  (24) 

CO  ( v)  + 02  (v  = 0)  -*■  CO  (v  - 1)  + 02(v  = 1).  (25) 


All  of  the  species  except  CO  are  assumed  to  be  entirely  in  their 
ground  vibrational  states,  so  that  the  reverse  processes  of  Eqs. 

(22)  - (25)  are  not  allowed.  The  CO-CO  and  C0-02  V-V  transfer 
probabilities  have  been  discussed  above  in  the  relaxation  calculations. 
The  CO-CS 2 / CO-CS,  and  CO-SO  collisions  were  treated  using  only  short 
range  forces  as  was  done  for  C0-02. 

For  the  V-T  processes,  the  general  event  is 

CO  ( v)  + M -*■  CO  (v  - 1)  + M,  (26) 

where  M = CO,  CS2,  CS , SO,  02,  He,  and  S.  All  of  these  V-T  tran- 
sition probabilities  have  been  computed  using  the  standard  SSH  theory. 
In  the  cases  CO-He  and  C0-02  the  computed  results  were  scaled  to 
match  experimental  data. 

There  is  one  further  difference  between  the  calculations  with 
chemical  kinetics  and  the  relaxation  calculations.  Inasmuch  as  the 
fluid  flow  evolves  at  constant  pressure,  and  therefore  the  gas 
kinetic  temperature  is  variable,  the  V-V  and  V-T  probabilities  are 
allowed  to  vary  with  temperature.  This  temperature  dependence  is 
described  parametrically  in  the  form: 
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r 


i 

I 

I 

I 

I 

: 


P(T)  = exp(a  + bT  + cT~l),  (27) 

where  the  parameters  a,  b,  and  c are  chosen  to  fit  the  computed 
probabilities  at  100  K,  300  K,  and  700  K. 


: 

I 

0 

.. 
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« 


I 
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2 . 3 Results 


The  kinetics  program  described  above  has  been  run  with  two  cases 
chosen  to  reproduce  as  closely  as  possible  the  experimental  conditions 
of  the  MDRL  CO  chemical  laser  devic  . . The  "standard  conditions" 
described  below  have  been  used  to  obtain  experimental  data  reported 
in  the  First  Quarterly  Progress  Report.^"  All  of  the  conditions 
necessary  for  these  two  runs  of  the  program  are  listed  in  Table  I, 
with  footnotes  as  additional  explanation. 

Turning  first  to  the  fluid  flow  and  chemical  kinetics,  Fig.  18 
shows  the  computed  stream  thickness,  temperature,  and  flow  velocity 
for  Case  I.  The  computed  quantities  at  x = 0.6  cm  are  of  interest 
since  the  cavity  optical  axis  is  located  at  that  position.  The.  only 
quantity  of  these  three  which  has  been  directly  measured  is  the  gas 
temperature.  An  epoxy-coated  thermocouple  positioned  at  x = 0.6  cm 
showed  a temperature  rise  with  the  addition  of  CS2  to  the  0,  02,  He 
stream  of  75  K,  while  the  calculated  rise  is  53  K.  This  point  will 
be  discussed  further  below;  the  values  of  stream  thickness  and  flow 
velocity  are  reasonable  and  agree  with  our  qualitative  experimental 
estimates . 

The  chemical  species  densities  are  shown  in  Fig.  19.  The  gas 
expansion  effects  are  included  in  these  densities;  the  He  density 
varies  from  1.806  x 1017/cm3  at  x = -1.2  cm  to  1.137  x 1017/cm’  at 
x = 10.0  cm.  In  these  cases  the  reactions  are  operated  oxygen  atom 
rich  on  the  basis  of  our  estimated  10%  02  dissociation,  and  the  CS2 
is  more  than  90%  reacted  at  x > 3.0  cm.  Consequently  the  CO  con- 
centration rises  to  a steady  state  value  of  2.45  x 1015/cm3, 
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corresponding  to  a partial  pressure  of  about  0.1  Torr.  The  difference 
between  the  CO  and  S densities  is  a measure  of  the  participation  of 
the  S + 0 2 reaction. 

Figures  20  through  24  show  the  population  density  and  optical 
gain  coefficients  for  Case  I.  Since  the  mixing  distance  is  taken 
to  be  0.5  cm,  the  initial  effect  of  the  reactions  is  to  build  up  an 
excited  distribution  of  CO  with  relatively  little  relaxation.  As 
Fig.  20  shows , the  CO  distribution  seen  by  the  optical  cavity  at 
x = 0.6  cm  (or  t - 37  psec)  will  be  close  to  the  CO*  distribution. 

By  x > 3 cm  (t  > 100  psec)  , the  CS 2 has  been  nearly  depleted,  and 
the  CO  distribution  begins  to  broaden  by  collisional  relaxation. 

Figures  21  and  22  show  the  resulting  maximum  P-  and  R-branch 
gains  versus  v with  time  (or  position)  as  a parameter.  Since  these 
are  zero-intensity  gain  coefficients,  the  only  comparison  to  be 
made  to  experiment  is  the  predicted  bands  versus  those  observed. ^ ^ 
Taking  the  optical  cavity  losses  to  be  a minimum  of  3%  in  the  single 
line  measurements,  bands  with  a value  of  a greater  than  3 x 10-4  cm-1 
(at  x = 0.6  cm)  should  be  seen.  From  Figs.  21  and  22  the  predicted 
bands  are  14  -*  13  down  to  6 5.  The  observed  bands  in  Table  II  are 

12  ->  11  down  to  7-6  (for  the  P-branch)  . This  suggests  that  the 
estimated  10%  dissociation  of  O2  is  too  high,  with  the  resulting  CO 
density  and  P-branch  gain  coefficients  too  large  in  magnitude.  A 
more  serious  difficulty  arises  with  the  R-branch  gain  coefficients 
of  Fig.  24.  The  maximum  a on  the  9 ->  8 band  is  predicti...  to  be 

only  a factor  of  4. 4/6. 4 lower  than  ap  at  x = 0.6  cm.  Similarly, 

_ 2 

gains  of  10  or  larger  are  predicted  for  ap  on  the  11  -*■  10  through 

39 


MCDOtVIVELL  DOUGLAS 


Fig.  20  CO*  population  distributions  for  Case  I.  Equivalent  times 
(cfter  CS2  injection)  and  positions  are  given. 


Fig.  21  Maximum  P-branch  gains  resulting  from  the  population  distributions 


Fig.  22  Maximum  P-branch  gains  versus  position  from  Fig.  20. 


Fig.  23  Maximum  R-branch  gains  from  Fig.  20. 


Table  II. 


Single  line  spectra  observed  in  MDRL  CO  chemical 
laser  with  conditions  of  Case  II  of  Table  I. 


Band,  branch 


y/  = line  observed 
0 = unable  to  get  oscillation 
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7 •+■  6 bands.  Our  experimental  observations  were  that  R-branch  lines 
could  be  seen  only  with  cold  CO  injection  (Case  II). ^ A scale 
change  in  the  CO  density  (and  hence  all  of  the  excited  CO  populations 
and  gain  coefficients)  will  not  resolve  this  discrepancy  since  the 
population  distribution  at  x = 0.6  cm  is  nearly  characteristic  of 
the  CO*  distribution  unmodified  by  relaxation. 

Figures  25  through  29  display  the  populations  and  gain  coef- 
ficients for  Case  II  which  is  identical  to  Case  I except  for  cold 
CO  injection  at  x = 0.0  cm.  (This  injection  leads  to  the  dis- 
continuities of  the  Up  and  aR  in  Figs.  28  and  29  because  of  changes 
in  CO*  densities  and  in  the  pressure-broadened  linewidths.)  The 
effects  of  the  cold  CO  injection  are  primarily  to  suppress  gain  on 
the  low  bands  at  long  times  (v  < 5 for  t i 250  sec),  and  to  enhance 
the  gains  on  higher  bands.  For  example,  at  x = C.6  cm,  a on  7 -+  6 
is  a factor  of  2.4  larger  with  cold  CO  injection  than  without. 

It  is  recommended  that  some  refinements  be  made  to  this  model 
of  the  CO  chemical  laser.  The  fluid  flow  description  should  be 
modified  to  include  viscous  effects  which  will  tend  to  broaden  and 
si  )w  the  gas  stream.  Additional  reaction  steps  may  be  added  to 
the  chemical  kinetics.  Perhaps  most  important,  we  suspect  the  CO* 
distribution  of  Reference  6 may  be  in  error  for  v < 6.  Other  work 
has  resulted  in  a distribution  which  has  substantial  formation  in 
the  states  below  v = 6.^  Changes  in  the  CO*  distribution  will  have 
profound  effects  on  the  stream  quantities  (through  heat  release  to 
translation  by  the  0 + CS  reaction) , and  on  the  optical  properties 
of  the  CO  medium.  The  CO*  distribution  may  resolve  the  discrepancy 
discussed  above  in  the  P-  and  R-branch  gain  coefficients  of  Case  I. 
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Fig.  25  CO*  population  distributions  for  Case  II  (with  CO  (300  K)  injection 


Fig.  28  Maximum  R-branch  gains  from  Fig.  25 


III. 
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